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I .  INTRODUCTION 

The  primary  objective  of  this  research  task  was  to  determine  the  displace¬ 
ment  range  of  laser  speckle  interferometry.  Interferometry  utilizes  light 
interference  phenomena  for  non-destructive,  sensitive  determinations  of  in¬ 
plane  displacements.  The  displacement  data  can  be  differentiated  to  obtain 
strains  which  can  in  turn  be  correlated  with  shearing  stresses  [1]. 

The  laboratory  setup  for  preparation  of  the  interferograms  utilized  an 
argon  laser  as  a  coherent  light  source  to  illuminate  the  specimen.  The  inter¬ 
ferograms  were  prepared  by  double  exposure  photography  techniques.  Analysis  of 
the  interferograms  employed  a  computer  aided  scanning  system  that  was  developed 
primarily  for  evaluation  of  speckle  interferometric  data. 

Conventional  techniques  of  measuring  surface  strain  and  displacement  data 
utilize  strain  gages,  dial  gages  and  various  other  electromechanical  sensing 
devices.  Although  such  methods  are  extremely  accurate,  information  is  limited 
to  very  small  regions  [2].  These  techniques  are  ideal  for  homogeneous  materials 
however,  nonconvent ional  applications  are  required  for  anisotropic  nonhomogen- 
eous  materials  such  as  composites.  Laser  speckle  interferometric  data  provides 
the  investigator  with  non-contact  full-field  information  over  an  entire  region 
of  interest. 

II.  OPTICAL  TEST  SETUP  AND  PROCEDURE 

The  experimental  test  arrangement  used  to  determine  the  interf erogram 
range  is  shown  in  Figures  1  and  2.  A  Spectra  Physics  argon  laser  adjusted  at 
0.50  watts  was  used  as  a  coherent  light  source  to  illuminate  the  specimen.  A 
double  exposure  photography  technique  using  AGFA  -  GEVAERT  Holotest  10E-56 
glass  film  plates  was  used  to  prepare  the  interferograms.  The  interferograms 
were  developed  with  conventional  photography  techniques. 

Figure  3  shows  the  specimen  setup  used  in  speckle  range  determination. 

The  translation  stage  was  positioned  on  a  vibration  free  table  by  means  of  a 
magnetic  base.  A  lens  positioner  assembly  was  mounted  on  the  translation 
table  with  2  polished  steel  rods  1.27  cm  in  diameter.  A  mounting  block 
extending  from  the  lens  positioner  contained  2  small  steel  rods  on  which  the 
rotation  stage  was  mounted.  The  specimen  was  attached  to  a  small  block  that  was 
bolted  on  the  rotation  table. 

Three  people  were  involved  in  preparing  a  complete  set  of  interferograms. 
One  person  rotated  the  specimen  between  each  exposure  by  means  of  a  micrometer 
control.  Another  person  on  the  optics  side  of  the  setup  shielded  the  lens 
between  each  exposure  so  as  not  to  destroy  the  interferograms.  A  third  person 
monitored  the  remote  shutter  control,  timer,  and  lights.  The  complete  inter- 
f erogram  preparation  process  was  done  in  the  dark.  Any  leakage  bf  white  light 
on  the  film  plate  would  destroy  the  interf erogram. 

A  double  exposure  photography  technique  was  utilized  to  make  the  inter¬ 
ferograms.  The  first  exposure  was  made  with  the  specimen  rotated  0.05°  from 
a  center  reference  point.  Next  the  camera  lens  is  covered  to  prevent  over¬ 
exposure  of  the  film  plates;  the  specimen  is  rotated  0.10°  and  a  second 
exposure  is  taken.  The  time  required  for  each  exposure  is  four  seconds. 
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Figure  1.  Experimental  test  setup  (sketch) 
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Figure  2.  Experimental  test  setup  (photo). 
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DATA  ANALYSIS 


A.  Computer  Scan  System 

The  experimental  displacement  data  was  collected  by  the  interfero¬ 
metric  analysis  system  shown  in  Figure  4.  This  system  utilizes  an  optical 
setup  and  computer  hardware  which  was  developed  primarily  for  applications  in 
speckle  interferometric  data  analysis.  The  inter ferograms  were  placed  in  a 
viewing  window  of  the  x-y  translation  table  and  reilluminated  by  a  HeNe  laser 
beam.  The  illumination  produces  a  fringe  pattern  on  a  diffuser  screen  which 
included  a  360  protractor  used  to  measure  the  displacement  angle  of  the 
fringes.  The  displacement  of  the  region  illuminated  is  inversely  proportional 
to  the  spacing  between  fringes.  The  direction  of  displacement  is  along  an 
axis  orthogonal  to  the  fringe  orientation.  A  typical  fringe  pattern  is  shown 
in  Figure  5. 

A  very  small  region  of  the  interf erogram  is  illuminated  and,  there¬ 
fore,  the  data  yields  the  displacement  at  a  point.  In  order  to  obtain  a 
complete  description  of  the  surface,  the  interf erograms  were  scanned  by  the 
x-y  translation  table  using  the  coordinate  system  shown  in  Figure  6.  The 
translation  stage  is  capable  of  0.0025  cm  incremental  movements  with  a 
0.00025  cm  repositioning  accuracy.  The  location  of  the  desired  point  on  the 
interferogram  was  controlled  by  translation  stage  stepping  motors  which  were 
in  turn  controlled  by  the  computer.  The  computer  program  used  to  control  the 
scan  pattern  and  calculate  the  displacements  appears  in  the  Appendix. 

B.  Displacement  Determination 

1.  Analytical  Considerations 

A  sketch  of  the  geometry  utilized  to  analytically  determine 
the  displacement  is  presented  in  Figure  7. 

Equidistant  rotation  of  length  T  about  a  center  point  produced 
an  isosceles  triangle  which  could  be  analyzed  with  the  law  of  sines  formulation 
as  given  by  Leithold  [3].  The  angle  of  rotation  was  determined  by  turning 
the  micrometer  7.9  microns  per  0.10  while  making  the  interf erograms .  With 
two  sides  and  an  angle  known,  the  displacement  can  now  be  determined. 

L  T 

SIN  a  =  SIN  8  (1) 

SIN  a 

so  L  =  SIN  8  T  (2) 

where  L  =  Displacement 

a  =  Angle  of  Rotation 

8  *=  180°  -  (90°  +  |) 

T  =  Length  of  axis  of  rotation 
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Figure  4.  Computer  analysis  system. 


Figure  7.  Geometrical  angular  rotation. 
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2.  Exper imental  Considerations 

A  computer  aided  scanning  system  was  utilized  to  analyze  the 
interf erograms  experimentally.  The  interferograms  were  placed  on  an  x-y 
translation  table  controlled  by  a  computer  program.  The  scanning  setup  used 
to  project  halos  of  the  interferograms  onto  a  diffuser  screen  is  shown  in 
Figure  8.  The  computer  controlled  translation  table  housing  the  interfero- 
gram  was  moved  in  a  raster  pattern  through  a  stationary  laser  beam.  Fringe 
spacing  and  the  angle  of  orientation  at  each  data  point  were  entered  into 
the  computer  which  in  turn  calculated  the  displacements  from: 


U 


0 


mi  f  s 
D 


(3) 


where 


U  =  Experimental  inplane  displacement 

D  =  Spacing  between  fringes 

m  =  Fringe  order 

f  =  Distance  from  interferogram  to 

analyzer  screen 
s  =  Film  scale  factor 

X  =  Wavelength  of  laser  light 


Ihe  complete  derivation  of  Equation  (3)  appears  in  Mullinix  [4].  The  vertical 
and  horizontal  components  of  displacements  were  also  determined  from: 


U 


V 


mXf  s 
D 


cos  0 


(4) 


U 


H 


mX  f  s 
D 


sin  0 


(5) 


The  wavelength  of  the  HeNe  laser  used  in  the  interferometric  analysis  is 
"128X.  The  fil  m  scale  factor  is  found  by  dividing  the  true  length  of  the 
specimen  by  the  image  length  on  the  interferogram.  The  computer  programs 
utilized  in  the  analysis  appear  in  the  Appendix. 

IV.  RESULTS 


The  primary  purpose  of  this  research  project  was  to  determine  the 
maximum  and  minimum  effective  range  of  laser  speckle  interferometry.  The 
maximum  range  determined  by  the  method  described  was  35149.0  t‘cm  with  an 
error  of  1.627.  from  the  analytical  value.  The  minimum  range  was  found  to  be 
''.85.8  ;  cm  with  an  error  of  3.037.  The  average  error,  or  total  coefficient 
of  variation  for  analytical  and  experimental  values  of  all  points  scanned 
(iables  1  through  3)  was  5.967..  This  error  could  be  attributed  to  dis¬ 
crepancies  in  angular  rotation  of  the  specimen,  inaccurate  measurements  taken 
lor  computer  scan  program,  and  reader  error. 

Horizontal  scan  plots  (Figures  9  through  12)  indicate  that  there  is  a 

direct  relationship  between  angle  of  rotation  and  amount  of  displacement. 

o  o 

I  ip, ure  9  shows  that,  as  the  angle  of  rotation  increases  from  0.10  to  0.20 
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TABLL:  1.  VKRTICAL  SCAN  DATA  SRD-T9  (0.10  ROTATION) 
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TABLE  2.  VERTICAL  SCAN  DATA  SRD-40  (0.15  ROTATION) 
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TABLE  3.  VERTICAL  SCAN  DATA  SRD-41  (0.20°  ROTATION) 
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at  a  fixed  Y  =  1.27  cm  value,  there  is  a  proportional  increase  in  the  dis¬ 
placements.  In  addition,  it  is  seen  that  the  displacements  are  symmetrical 
about  the  center  line  of  the  specimen,  which  is  located  where  the  fringe  angle 
is  zero.  Data  for  these  plots  are  displayed  in  Tables  4  through  9. 

V.  SUMMARY  AND  CONCLUSIONS 

The  laser  speckle  interferometric  range  was  determined  by  applying  an 
optical  process  that  utilizes  a  coherent  laser  light  source  for  precise 
determination  of  linear  inplane  displacements.  The  analysis  of  the  interfero- 
grams  was  accomplished  by  a  computer  aided  scanning  system  designed  primarily 
for  evaluating  interferometric  data.  This  nondestructive,  highly  sensitive 
method  of  analysis  is  extremely  accurate  and  provides  the  investigator  with 
full-field  displacement  data.  A  full-field  reconstructed  optical  Fourier 
transform  image  of  the  displacement  region  can  also  be  produced  by  this 
cechn ique . 


For  this  research  task,  a  composite  test  specimen  was  scanned  both 
vertically  and  horizontally  at  four  specified  regions.  The  maximum  and 
minimum  displacement  was  determined  to  be  35149. Oticm  and  685.8 ucm ,  respectively. 
Horizontal  scan  evaluation  indicated  that  the  rotational  angle  and  the  amount 
of  displacement  increase  proportionally.  This  research  exercise  demonstrates 
the  capability  to  capture  deformation  under  load  of  rotating  objects  such  as 
helicopter  blades  and  high  speed  turbine  vanes.  A  pulsed  high  energy  coherent 
light  source  and  high  resolution  multiple  exposure  film  would  be  required  to 
capture  the  interferometric  data.  The  inter f erogram  could  be  investigated 
wotli  qualitatively  for  anomalies  in  the  structures  as  well  as  quantitatively 
to  provide  deformation  under  load  profiles  for  engineering  design  applications. 

Attempts  have  been  made  by  Swinson  [5]  to  increase  the  effective  range 
of  laser  speckle  interferometry  coupled  with  accoustical  interferometry  through 
sandwiching  techniques.  The  technique  introduces  a  rigid  body  shift  into  the 
deformed  state  of  a  test  object  to  increase  the  number  of  readable  fringes; 
after  which,  the  rigid  body  effects  are  subtracted  out  of  the  analysis. 

There  is  a  definite  need  for  rapid  advancement  in  the  state-of-the-art  of  laser 
speckle  interferometry  to  include  the  most  advanced  methods  of  computer 
analysis  and  optical  image  processing. 
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-1.524  [  6943  4416  -5357  -50.5  15200  14190  -5447  -21. 


TABLE  9.  HORIZONTAL  SCAN  DATA  SRD-41  (0.20°  ROTATION) 


\  Y 
x\ 

(cm)\ 

6.35  cm 

U 

(pcm) 

uH 

(pern) 

UV 

(pc  m) 

0 

1.524 

25563 

24957 

5533 

12.5 

1.270 

23433 

23040 

4270 

10.5 

1.016 

21630 

21364 

3384 

9.0 

0.762 

20085 

19935 

2448 

7.0 

0.508 

18746 

1867  5 

1634 

! 

5.0 

0.254 

18746 

18743 

327 

1.0 

0.000 

18746 

18743 

-  327 

-  1.0 

— 

-0.254 

19393 

19356 

-1184 

-  3.5 

-0.508 

19393 

19303 

-1859 

-  5.5 

-0.762 

21630 

21512 

-2261 

-  6.0 

-1.016 

22495 

22218 

-3519 

-  9.0 

-1.270 

23433 

22962 

-4672 

-11.5 

-1.524 

24452 

23383 

-7149 

-17.0 

NOTE:  Data  at  8.89  cm  out  of  range. 
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APPENDIX 
COMPUTER  CODES 
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Program  Name: 
Program  Function: 


SCAN 


Input : 


Output : 


To  evaluate  interferograms,  compute  horizontal  and 
vertical  components  of  displacement  as  well  as  total 
displacement . 

1)  Stage  to  be  advance  -  x-direction  or  y-direction 

2)  Stage  Increment  -  Distance  moved  between  scan  points 

3)  Film  scale  factor  -  Ratio  of  real  object  to  image  on 
interf  erogram 

4)  Distance  to  diffuser  screen  -  Distance  from  interf erogram 
to  diffuser  screen 

5)  Displacement  -  Distance  measured  manually  between  fringes 

6)  Angle  -  Angle  of  orientation  of  fringes 

1)  -  Total  Distance  from  Initial  Scan  Point 

2)  -  Total  Displacement 

3)  -  Horizontal  Component  of  Displacement 

4)  -  Vertical  Component  of  Displacement 


i 


COMMON  DD  (2,200) 

WRITE(5,1) 

1  FORMAT  ('MANUAL  YOUNG  FRINGE  ANALYZER') 

WRITE(5,2) 

2  FORMAT  ('STAGE  TO  BE  ADVANCED?  0=X,  1=Y)') 

READ (5, 3)  IA 

3  FORMAT (II) 

WRITE(5,4) 

4  FORMAT ('STAGE  INCREMENT?') 

READ (5, 5)  IC 

5  FORMAT (14) 

WRITE(5,6) 

6  FORMAT( 'FILM  SCALE  FACTOR?') 

READ(5,7)  SF 

7  FORMAT (F 10 . 0) 

WRITE(5,17) 

17  FORMATC DISTANCE  FROM  INTERFEROGRAM  TO  DIFFUSER  SCREEN?') 
READ(5,7)  X 
ICC=0 

8  ICC=ICC+1 
ICP=ICC-1 
WRITE (5, 9) 

9  KORMAT( 'DISPLACEMENT? ') 

READ (5, 7)  D 

WRITE(5, 10) 

10  FORMAT ('ANGLE?') 

READ (5,7)  A 
A=A*3. 14159/180. 

IF(D.EQ.O)  GOTO  13 
IE(D.LT.O)  GOTO  15 
U=(SF*.Q000249133*X)/D 
Ul=U*COS (A) 

U2=U*SIN(A) 

H=ELOAT ( ( ICC-1) *IC) * . 001 
WRITE(5,11)  H,U ,U1 ,U2 

11  FORMAT ( 311  H=  , FI 0 . 3 , 5X ,  2HU= ,  F12 . 6) 

13HU1 = , F12 . 6, 5X, 3HU2= ,F12 . 6) 

12  DD ( 1 , 1 CC)=U 1 
DD(2 , I CC)=U2 
GOTO  14 

13  DD( 1,1 CC)=0 . 

1)D(2 , 1CC)=0. 

14  IF(lA.EQ.l)  CALL  YADV(IC,5) 

IF(IA.EQ.O)  CALL  XADV(IC,5) 

GOTO  8 

15  ICC=(1-ICC)YIC 
IF(IA.EQ.l)  CALL  YADV(ICC,5) 

IF(IA.EQ.O)  CALL  YADV(ICC,5) 

WRITE(5 , 16) 

16  FORMAT ( 'ANALYSIS  ENDED') 

!>  1  OP 
END 
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SCAN  PROGRAM  (CONTINUED) 


SUBROUTINE  YADV(IS.IR) 

IS=NO.  STEPS  (+**FWD ,  -=REV) 
IR=ADVANCE  RATE  OF  STAGE 
X=0 

IF(IS.GT.O)  GOTO  3 
IF=IABS(IS) 

DO  2  1=1, IP, 1 

CALL  IPOKE("167772, ”020000) 

DO  7  K=l,  IR,1 

YSIN(X)CALL  I POKE  ("167772, "000000) 

DO  1  J=1,IR,1 

Y=SIN(X) 

CONTINUE 

GOTO  6 

CONTINUE 

DO  5  11=1, IS, 1 

CALL  IPOKE("167772, "010000) 

DO  8  KK=1,IR,1 
Y=SIN(X) 

CALL  IPOKE("167772, "000000) 

DO  4  JJ=1 , IR, 1 
Y=SIN(X) 

CONTINUE 

CONTINUE 

RETURN 

END 
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SCAN  PROGRAM  (CONTINUED) 


SUBROUTINE  XADV(IS.IR) 
IS=NO.  STEPS  (+=FWD,-=REV) 
IR=ADVANCE  RATE  OF  STAGE 
X=0 

IF(IS.GT.O)  GOTO  3 
IP=IABS(IS) 

DO  2  1=1, IP, 1 

CALL  IPOKE("167772, "100000) 

DO  7  K=1 , IR, 1 

Y=SIN(X) 

CALL  IPOKE("167772, ”000000) 

DO  1  J=1 , IR, 1 

Y=SIN(X) 

CONTINUE 

GOTO  6 

CONTINUE 

DO  5  11=1, IS, 1 

CALL  IPOKE("167772, "040000) 

DO  8  KK=1 , IR , 1 

Y=SIN(X) 

CALL  IPOKE("167772, "000000) 

DO  4  J J=1 , IR, 1 

Y=SIN(X) 

CONTINUE 

CONTINUE 

RETURN 

END 
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Program  Name:  Move 

Program  Function:  To  move  translation  stage  and  interf erogram  to  initial 
scan  position 

Input:  Input  X,Y  -Four  digit  value  for  positive  or  negative 

movement  in  X  and  Y  directions 

Output:  Motion  of  translation  stages 

0.0025cm  Incremental  Movements 
0.00025cm  Reposition  Accuracy 

3  WRITE(5,1) 

1  FORMATC INPUT  X,Y-214') 

READ(5,2)  IX, IY 

2  FORMAT (2 14) 

IF(IX.EQ.O. AND . I Y . EQ . 0)  GOTO  4 
CALL  XADV(IX,5) 

CALL  YADV(IY,5) 

GOTO  3 

4  CONTINUE 
STOP 
END 

SUBROUTINE  YADV(IS.IR) 

C  IS=N0 .  STEPS  (+=FWD,-=REV) 

IR=ADVANCE  RATE  OF  STAGE 
X=0. 

IF(IS.GT.O)  GOTO  3 
IP=IABS (IS) 

DO  2  1=1, IP, 1 

CALL  IPOKE("167772, "020000) 

DO  7  K=1 ,IR,1 

7  Y=SIN(X) 

CALL  IPOKE("167772, "000000) 

DO  1  J=1,IR,1 

1  Y=SIN(X) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  II-1,IS,1 

CALL  IPOKE("167772, "010000) 

DO  8  KK=i , IR, 1 

8  Y=SIN(X) 

CALL  IPOKE( "167772, "000000) 

DO  4  JJ=1,IR,1 

4  Y-SIN(X) 

5  CONTINUE 

6  CONTINUE 
RETURN 
END 

SUBROUTINE  XADV(IS,IR) 

C  IS=NO.  STEPS  (4-FWD.--REV) 


37 


MOVE  PROGRAM  (CONTINUED) 


C  IR=ADVANCE  RATE  OF  STAGE 
X=0 

IF(IS.GT.O)  GOTO  3 
IP=IABS(IS) 

DO  2  1=1, IP, 1 

CALL  IPOKE("167772, ”100000) 
DO  7  K=1 , IR, 1 

7  Y=SIN(X) 

CALL  IPOKE("167772, ”000000) 
DO  1  J=1 , IR, 1 

1  Y=SIN(X) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  11=1, IS, 1 

CALL  1P0KE("167772, ”040000) 

DO  8  KK=l , IR, 1 

8  Y=SIN(X) 

CALL  IPOKE("167772, ”000000) 
DO  4  JJ=1,IR,1 

4  Y=S IN (X) 

5  CONTINUE 

6  CONTINUE 
RETURN 
END 
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